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ABSTRACT 

We present the Final Analysis of the European Large Area ISO Survey (ELAIS) 15 
/zm observations, carried out with the ISOCAM instrument on board the Infrared 
Space Observatory (ISO). 

The data reduction method, known as LARI method, is based on a mathematical 
model of the detector's behaviour and was specifically designed for the detection of 
faint sources in ISO-CAM/PHOT data. The method is fully interactive and leads to 
very reliable and complete source lists. 

The resulting catalogue includes 1923 sources detected with S/N > 5 in the 
0.5 - 100 mJy flux range and over an area of 10.85 deg^ split into four fields, making 
it the largest non-serendipitous extragalactic source catalogue obtained to date from 
ISO data. 

This paper presents the concepts underlying the data reduction method together 
with its latest enhancements. The data reduction process, the production and basic 
properties of the resulting catalogue are then discussed. 

Key words: infrared: galaxies - galaxies: formation, evolution, active, starburst - 
cosmology: observations - methods: data analysis - catalogues. 



1 INTRODUCTION 

The Infrared As t ronomical Sat ellite (IRAS, 
iNeueebauer et al.l Il984l: ISoifer et al.l Il98'it was extremely 

* e-mail: m.vaccari@imperial.ac.uk 



successful in characterizing for the first time the global 
properties of the mid- and far-infrared sky, carrying out 
an all-sky survey at wavelengths of 12, 25, 60 and 100 /im 
and leading to discoveries such as those of Luminous, Ul- 
traluminous and Hyperluminous Infrared Galaxies (LIRGs, 
ULIRG and HLIRGs, respectively), a substantial popula- 
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tion of evolving starbursts and the detection of large-sca le 
structure in the galaxy distribution dSaunders et al.lll99J) . 

Unfortunately, the IRAS view was typically limited to 
the very local Universe {z ^ 0.2), thus hampering statistical 
studies of infrared-luminous galaxies at cosmological red- 
shifts. Only few sources were detected by IRAS at higher 
redshifts, typically ULIRGs mag nified by gravitational 
lenses, like F10214-(-4724 [z = 2.28, 'Rowan-R obinson et alJ 
llQQlT) . In particular only about 1000 galaxies were de- 
tected all over the sky in IRAS 12 /xm band. Infrared 
sourc e counts base d on I RAS data iRowan-Robinson et alJ 
Il984l: ISoifcr ct al.^ Il984^ showed some marginally signif- 
icant excess of faint sources with respect to no evolu- 
tion models ([H acking et al] Il987l : iFrances jhmiet.^ 
iLonsdale et alJ l^gOtlGTCgori^h^^^ 119951: iBertin et alJ 
Il997r . but not enough statistics and dynamic range in flux to 
discriminate between evolutionary scenarios were available. 

Although conceived as an ob servatory-type miss ion, the 
Infrared Space Observatory (ISO. lKessler et alJiggsT) was in 
many ways the natural successor to IRAS, bringing a gain of 
a factor ~ 1000 in sensitivity and ~ 10 in angular resolution 
in the mid-infrared. A substantial amount of ISO observing 
time was therefore devoted to field surveys aimed at detect- 
ing faint infrared galaxies down to cosmological distances. 
Such surveys were conceived as complementary in flux depth 
and areal coverage, allowing a systematic investigation of 
the extragalactic sky down to so far unattainable flux densi- 
ties at both mid an d far infrared wavelengths, w hose results 
are summarized bv lGenzel a nd Cesa rskvl i2000t) . In partic- 
ular, extr agalactic 15 /xm s ource c ounts determi ned with 
IS OGAM jElbaz et al.ll999l:rGruppioni et al.l2002D have re- 
vealed a significant departure from Euclidean slope within 
the 1-5 mjy flux range, which has been interpreted as evi- 
dence for a strongly evolving population of starburst galax- 
ies. 

The Europea n Large Area ISO Survey (ELAIS, 
lOliver et al] I2OO0I: iRowan- Robinson et al] l2003f) was the 

most ambitious non-serendipitous survey and the largest 
Open Time project carried out with ISO, aimed at bridg- 
ing the flux gap between IRAS all-sky survey and ISO 
deeper surveys. ELAIS observations mapped areas of about 
12 deg^ at 15 and 90 l ira and smaller areas at 7 and 175 
/xm with the IS OGAM (|Gesarskv et al.lll996l . 7 and 15 ^m) 
and ISOPHOT (Lemk e et al.lll996l . 90 and 175 ^m) cam- 
eras. Most importantly, ELAIS 15 /xm observations are the 
only ones allowing to sample the 1-5 mJy flux range, where 
most of the source evolution appears to take place. 

Since the project approval, the ELAIS consortium, 
grown in time to a total of 76 collaborators from 30 Eu- 
ropean institutes, has undertaken an extensive program of 
ground-based optical and near-infrared imaging and spec- 
troscopy. Thanks to such an extensive multi- wavelength cov- 
erage, the ELAIS flelds have now become among the best 
studied sky areas of their size, and natural targets of on- 
going or planned large-area surveys with the most powerful 
ground- and space-based facilities. Further details on ELAIS 
multi- wavelength observations and catalogues are presented 
in lRowan- Rob inson et al. (2003). After the loss of the WIRE 
satellite, notwithstanding the observations at several in- 
frared wavelengths soon to come from Spitzer and later from 
SOFIA and Herschel, ISO observations will remain a valu- 
able database for many years to come. In particular, until 



the advent of JWST, ELAIS 15 /xm observations will provide 
a complementary view on three areas (SI, Nl and N2) which 
will be covered at different wavelengths as part of the Spitzer 
Wide-Area Extragalactic Survey (SWIRE, . Lonsdale et alJ 
I2OO3). Thus the need of reducing such data with the utter- 
most care and provide the community with an agreed-upon 
legacy from the ELAIS project. 

This paper presents the Final Analysis of ELAIS 15 
^m observations, and is structured as follows. In Section |21 
a brief description of the most relevant aspects of ELAIS 
15 /im dataset is given. Section |21 describes the data reduc- 
tion method and its improvements. In Section 2] the tech- 
nique employed for flux determination and its results are 
presented. Section |3 details the results of the simulations 
that were carried out in order to assess the performance 
of the data reduction method and thus the quality of the 
resulting catalogue. In Sections [7| and |H| respectively, esti- 
mates of the achieved astrometric and photometric accuracy 
are given. Section |S| summarizes the identiflcation of 15 /xm 
sources in optical and near-infrared images, while Section |5| 
describes the procedure adopted to establish the catalogue 
photometric calibration. Finally, Section [TUl describes gives 
a basic description of the catalogue contents. 



2 THE ELAIS 15 DATASET 

The ELAIS 15 /xm main^ dataset is made up of 28 rasters 
(ISO basic imaging observations), each covering an area of 
about 40' X 40', divided into 4 fields, one (SI) in the southern 
hemisphere and three (Nl, N2 and N3) in the northern one. 
Small superpositions at the boundaries and a limited degree 
of redundancy on portions of the fields give a total covered 
area of 10.85 deg^. 

The fields were selected on the basis of their high eclip- 
tic latitude (|/3| > 40°, to reduce the impact of Zodiacal 
dust emission), low cirrus emission (/looMm < 1-5 MJy/sr) 
and absence of any bright {Si2nm > 0.6 Jy) IRAS 12 /xm 
source. In Figure the location on the sky of the survey 
fields is showU j overlaid on cirrus m aps (GOBE normalized 
IRAS maps of ISchlegel et alJ Jl998fl '). Nearby IRAS sources 
with 12 /im fiuxes brighter than 0.6 Jy are also plotted. The 
overall sky coverage, highlighting the position and redun- 
dancy of single rasters, is illustrated in Figure |21 

ELAIS 15 pim observations were carried out operating 
the ISOGAM instrument in raster mode using the LW3 fil- 
ter. In this observing mode, the ISOGAM 32 x 32 pixel LW 
detector was stepped across the sky in a grid pattern, with 
about half detector width steps in one direction and the 
whole detector width steps in the other. Thus, reliability 
was improved as each sky position (apart from those at the 
boundaries of the raster region) was observed twice in suc- 
cessive pointings and overheads were reduced because each 
raster covered a large area (~ 40' x 40'). At each raster 
pointing (i.e. grid position) the detector was read out several 

^ Note that smaller sky regions observed with an higher redun- 
dancy as part of the ELAIS project, such a s the S2 field (who se 
data reduction and analysis is described bv lPozzi et al 
the Xl,...,6 fields and the ultra-deep portion of the Nl field, which 
was observed ten times, are not considered in this work. See 
lOliver et alj i200(i) for further details on these smaller fields. 
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Figure 1. Sky Locations of ELAIS 15 fim Fi elds. From top left to bottom right: Nl, N2, N3 and SI areas. Greyscales indicate COBE 
normalized IRAS 100 /im intensity maps from lSchlegel et al.l il998ft . IRAS sources with 12 fim fluxes brighter than 0.6 Jy are also drawn 
as circles with radii proportional to their fluxes 



(typically 10) times, to increase the redundancy in order to 
be able to identify cosmic rays impacts and distinguish their 
severe effects on the electronics from real sources. Further- 
more, on the raster first pointing, 80 readouts were carried 
out to allow the detector to approach stabilization. Table Q 
describes the observation parameters for the LW3 observa- 
tions, while Table |21 hsts the fields and related ISO Data 
Archive filenames making up the dataset. 



3 DATA REDUCTION 

Reduction of data obtained with ISO instrumentation has 
always proved very difficult for a number of reasons. As far 
as ISOCAM observations carried out using its Long Wave- 



length (LW) detector are concerned, the two most important 
instrumental phenomena one has to deal with are the qual- 
itatively very different effects produced on the detector's 
electronics by the frequent and severe cosmic ray impacts, 
which have long been known and referred to as glitches, 
and its sizable transient behaviour after changes in the inci- 
dent photon flux, which we will hereafter simply refer to as 
transients. In both cases, the cryogenic operational tem- 
peratures of the detector caused it to very slowly respond 
after these events. Lack of an accurate modeling of these 
effects can thus lead to spurious detections or errors in flux 
determination. In ELAIS 15 data reduction, the impact 
of these effects is increased by the instrumental parameters 
which were chosen in order to maximize the survey area. 
The short integration time (2.1 s x 10 frames/pointing), the 
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Figure 2. ELAIS 15 /im sky coverage. Prom top left to bottom right: Nl, N2, N3 and SI fields. Note that different fields are represented 
with a slightly different image scale. 



Table 1. ELAIS ISOCAM LW3 Observation Parameters 



Parameter Value 



Band ^eff li.?) fim 

Band FWHM Range 12.0-18.0 /xm 

Detector gain 2f 

Integration time 2 s 

Number of exposures per pointing 10 

Number of stabilization exposures 80 

Pixel field of view 6" 

Number of pixels 32 X 32 

Number of horizontal and vertical steps 28 , 14 

Number of rasters (including repetitions) 28 

Horizontal and vertical step sizes 90" , 180" 

Total area 10.85 deg^ 



t Except in "test" raster N2_R_A where gain was 1 

large raster step (half the detector size along one axis and 
the whole detector size along the other one) and the large 
pixel size (6 arcsec) all contribute to reduce the redundancy 
and to increase the undersampling. Low redundancy and 
high undersampling, in turn, increase difficulties in distin- 
guishing sources from strong glitches (low reliability) and in 
correctly determining source fluxes (low photometric accu- 
racy). 

Roughly speaJting, glitches can be divided into three 
categories according to the way they shape the detector's 



output signal, their decay time and influence on the pixel 
responsivity: glitches belonging to these different classes are 
respectively dubbed common glitches, faders and dip- 
pers. Slow decreases of the signal following cosmic ray 
impacts are called faders, while prominent reductions of 
the pixel responsivity very slowly recovering afterwards are 
called dippers. These two effects are believed to be asso- 
ciated with proton or a particle impacts on the detector, 
and have a fairly long lasting impact on pixel responsivities. 
Conversely, the much more frequent impacts of cosmic ray 
electrons produce common glitches characterized by a rela- 
tively fast decay time, lasting only a few readouts. Therefore, 
the number of frames affected by a single fader or dipper is 
much higher than in the case of a common glitch, the pixel 
responsivity taking from tens to hundreds of seconds to re- 
cover completely. However, common glitches are much more 
frequent than faders and dippers and, may all the same ham- 
per the quality of data reduction. Thus, all kinds of glitches, 
if not correctly removed (or, more properly, corrected for), 
can lead to spurious detections, or unreliability, and to the 
loss of genuine sources, or incompleteness. Ou the other 
hand, transients all follow the same pattern, due to the fact 
that they arise from the non negligible time it takes for the 
output signal to reach the stabilization value after a change 
in the incident photon flux has taken place. The measured 
signal is thus always lower than the true one. Failing to 
model this time effect in data reduction can lead to a sys- 
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Table 2. ELAIS 15 /xm Fields. Field Name, Raster Name, ISO 
Data Archive (IDA) Official Filename and Raster Centre Coordi- 
nates 



Field 


Raster 


IDA 


RA (J2000) 


Dec (J 2000) 


SI 


S1.1 


23200251 


00 


30 


25.40 


-42 57 00.3 




S1.2 


23200353 


00 


31 


08.20 


-43 36 14.1 




S1.3 


41300955 


00 


31 


51.90 


-44 15 27.0 




S1.4t 


23300257t 


00 


33 


59.40 


-42 49 03.1 




S1.5.A 


23300459 


00 


34 


44.40 


-43 28 12.0 




S1.5.B 


77500207 


00 


34 


44.40 


-43 28 12.0 




S1.5.C 


78502406 


00 


34 


44.40 


-43 28 12.0 




S1.6t 


41001161t 


00 


35 


30.40 


-44 07 19.8 




81.7 


40800663 


00 


37 


32.50 


-42 40 41.2 




S1.8 


40800765 


00 


38 


19.60 


-43 19 44,5 




S1.9 


41001867 


00 


39 


07.80 


-43 58 46.6 


Nl 


Nl_l 


30200101 


16 


15 


01.00 


+54 20 41.0 




N1.2.A§ 


30400103§ 


ID 


1 Q 

io 




+54 59 35.9 




Nl.2^ 


67200103 


16 


13 


57.10 


+54 59 35.9 




N1.3 


30500105 


16 


10 


34.90 


+54 11 12.7 




N1.4 


30600107 


16 


09 


27.00 


+54 49 58.7 




N1.5 


31000109 


16 


06 


10.80 


+54 01 08.0 




N1.6 


30900111 


16 


04 


59.00 


+54 39 44.3 


N2 


N2.1 


50200119 


16 


32 


59.80 


+41 13 33.2 




N2_2 


51100131 


16 


34 


44.50 


+40 38 45.0 




N2.3 


50000723 


16 


36 


05.50 


+41 33 11.8 




N2_4 


50200225 


16 


37 


48.90 


+40 58 13.1 




N2.5 


50100727 


16 


39 


13.80 


+41 52 31.6 




N2_6 


50200429 


16 


40 


55.50 


+41 17 22.7 




N2.R.A 


11600721 


16 


35 


45.00 


+41 06 00.0 




N2.R.B 


77900101 


16 


35 


45.00 


+41 06 00.0 


N3 


N3.3 


42500237 


14 


29 


38.30 


+33 24 49.6 




N3.5 


43800341 


14 


32 


38.20 


+33 11 10.3 



t Originally flagged as "Telemetry Drops" in observation logs 
X Originally flagged as "Unknown Quality" in observation logs 
§ Originally flagged as " Aborted" in observation logs 



tematic underestimation of source fluxes. For these reasons, 
the data cleaning and modeling is an extremely delicate pro- 
cess requiring great care in order to produce highly reliable 
sky maps and source lists. 

While it was variously demonstrated that it is possible, 
at least to a certain extent, to describe the detector's be- 
haviour adopting some physical model, the large number of 
readouts involved in raster observations and the peculiar na- 
ture and strength of noise patterns also require efficient and 
robust algorithms to be developed so as to make the actual 
data reduction undertaking feasible in a nearly-automatic 
way. 

A number of data reduction methods has thus been de- 
veloped an d tested, mostly on deep fields (e.g. the PRETI 



Starck et al 



Desert et al 



199£ ) and the Triple Beam Switch 



1 19991) V Unfortunately, such meth- 



method by 
method by 

ods did not prove completely reliable on shallower fields, 
leading to a high number of false detections. Besides, these 
methods suffered from the lack of an efficient way to in- 
teractively check the quality of the data reduction when 
needed. The Prelim inary Analysis of ELAIS 15 ^tm data 
jSerieant et al.ll2000l) was thus carried out adopting a more 
traditional approach involving the corroboration of auto- 
mated detections through visual inspection by different ob- 
servers. 

The data reduction described in this paper was car- 



ried out using the LARI method jLari et all 1200 ll. l20o3) . 

a new technique developed to overcome these difficulties 
and provide a robust interactive technique for the reduction 
and analysis of ISOCAM and ISOPHOT data, particularly 
suited for the detection of faint sources and thus for the full 
exploitation of their scientific potential. The method was 
variously refined, and significantly better resul ts are now 
obtained, with respect to the technique used in lLari et alJ 
(2001) for the reduction of the SI field, so that a thorough 
rc-rcduction of all ELAIS fields seemed appropriate and was 
thus carried out. As before, data reduction is carried out 
within an IDL environment using mostly p urpose-built rou- 
tines, exploiting CAM Interactive Analysis dOtt et alJl200ll . 
CIA,) software for basic operations only. 



3.1 The Model 

The LARI method describes the sequence of readouts, or 
time history, of each pixel of ISOCAM LW3 detector in 
terms of a mathematical model for the charge release to- 
wards the contacts. Such a model is based on the assump- 
tion of the existence, in each pixel, of two charge reservoirs, 
a short-lived one Qt, also known as breve, and a long-lived 
one Qi, also known as lunga, evolving independently with 
a different time constant and fed by both the photon flux 
and the cosmic rays. The model is fully charge-conservative, 
i.e. no decay of accumulated charges is considered, except 
towards the contacts, and thus the observed signal 5* is re- 
lated to the incident photon flux / and to the accumulated 
charges Qt and Qi by the 

dQtot r dQ; 

dt ' 



dt ' dt dt ' 
where the evolution of these two quantities is governed by 
the same differential equation, albeit with a different effi- 
ciency Ei and time constant Oi 

dQ, 
dt 
so that 



at ( 



where i — b,l 



(2) 



S = [1 - et - ei) I + at 



(3) 



The values of the parameters d and a; depend on the de- 
tector's physics, and are thus assumed to depend on the ex- 
posure time of a given observation and on the stabilization 
background level of a given pixel, according to the relation 



t / S 



(4) 



where Uifi is the value of ai relative to a reference exposure 
time to and average signal level So. In practice, an additive 
offset signal attributable to thermal dark current (a compo- 
nent which is otherwise removed in standard dark current 
subtraction), is actually added to both S and I in the equa- 
tion above when it is estimated to be important, i.e. in the 
rare cases when the deepest dippers' depth otherwise ex- 
ceeds 10 % of the stabilization background level. 

3.2 Pipeline 

Data reduction begins with a preliminary pipeline incor- 
porating all necessary steps in order to prepare the data for 
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the temporal fitting procedure wliicli is tlie critical step of 
the LARI method. 

Raw data downloaded from the ISO Data Archive 
are first imported into the IDL raster structure contain- 
ing all observational information using CIA routines. Like- 
wise, dark current subtraction and conversion from ADU 
to ADU/gain/s are carried out using CIA. Then a dedi- 
cated IDL structure called liscio is built in order to con- 
tain not only all raster information but also all ancillary 
arrays needed in order to carry out the following reduction. 

At this stage potential glitches are identified through a 
two-step median filtering process of each pixel time history, 
or deglitching. Then a separate routine determined the sta- 
bilization background (or global background, as opposed to 
the local background described below) and the aforemen- 
tioned offset signal. Such routine was carefully devised to 
filter out in so far as possible long-term effects such as sta- 
bilization, faders and dippers and thus provide a reliable 
estimate of these two parameters the fitting procedure is 
particularly sensitive to. The same routine also identifies 
the few potential bright sources which at times the fit failed 
to recognize by itself, which are then interactively checked 
by eye to assess their reliability. 

3.3 Fitting 

The signal as a function of time is finally processed inde- 
pendently for each pixel. The fitting procedure models the 
transients attributable to changes in incident fiux and the 
features on both short and long time scales produced by 
cosmic ray impacts on the time history, modeling glitches as 
discontinuities in the charge release. As seen in Section f3. II 
the same values for the d and Ui parameters are used for all 
pixels, apart from the scaling of the ai according to the ex- 
posure time and the stabilization background level, leaving 
as free parameters only the charges at the beginning of the 
observation and at the "peaks" of glitches. 

Glitches from nearby pixels are also considered when 
their height is substantially (i.e. 20 times) higher than 
the chosen threshold. The fitting algorithm starts with the 
strongest potential glitches identified in pipeline deglitch- 
ing, assumes discontinuities at these positions and tries to 
find a fit to the time history that satisfies the model as- 
sumed to describe the solid-state physics of the detector. 
If no acceptable fit is found, the next fainter glitch is con- 
sidered as a potential discontinuity, and so on. Iteration of 
the fitting procedure is interrupted when either a satisfac- 
tory (typically 0.2 ADU/gain/s) data-model rms deviation 
is achieved or the maximum number of allowed iterations is 
reached. 

At this stage the code estimates several quantities 
needed to build the sky maps on which source extraction will 
then be performed. All of these quantities are "recovered" 
from glitches, i.e. already take into account discontinuities 
in charge release assumed to describe glitches during fit, and 
their list includes: 

• the charges stored into the breve and lunga reservoirs 
at each readout. 

• the local background, i.e. the signal to be expected on 
the basis of the previously accumulated charges if only the 
stabilization background were hitting the detector. 



• the model signal, produced by the incident flux com- 
ing from both the stabilization background and detected 
sources. 

• the "unreconstructed" signal, i.e. the detected fiux 
recovered from glitches but not from transients. This is com- 
puted as the difference between the measured signal and the 
local background. 

• the "reconstructed" signal, i.e. the detected fiux re- 
covered not only from glitches but also from transients as- 
sociated with changes in incident fiux. 

For the sake of clarity, one must emphasize here the differ- 
ences between the two kinds of signals (and the correspond- 
ing fiuxes and sky maps they will finally turn into) defined 
above, namely unreconstructed and reconstructed sig- 
nals. Both quantities take into account the effects of the sta- 
bilization background and glitches on the detector, and their 
values are thus expected to be negligible if a source is not 
illuminating the pixel at the given pointing. The difference 
between the two quantities only appears when an additional 
signal the code is not able to model otherwise (i.e. as the 
effect of a glitch) is detected and attributed to a source of 
a given flux. The code then models the transients expected 
from such an additional flux and "reconstructs" the signal 
one would detect if they did not affect the detector, i.e. if 
its response were instantaneous. In other words, unrecon- 
structed signals do not take into account the effects of this 
modeling of changes in incident flux, thus representing the 
effective charge collected during the exposure, whereas re- 
constructed signals recover the charge "loss" due to the slow 
detector response. Therefore, the former are systematically 
lower than the latter. Actually, simulations show that the 
code is not actually able to properly model transients be- 
low a certain intensity threshold, thus suggesting the use of 
unreconstructed signals only to carry out further processing. 

Figure 121 shows how a successful flt is able to describe 
cosmic-ray-induced violent changes in the signal level and 
thus recover useful information (speciflcally, source fluxes) 
from the pixel time history. Panels a) and b) show two ex- 
amples of how glitches (a fader and a dipper, respectively) 
are described as discontinuities in the signal level slowly re- 
covering towards the stabilization background, while panels 
c) and d) show how sources are detected even on the top of 
strong glitches. The solid line represents the observed data, 
the short-dashed line the best-fit model and the long-dashed 
line the detected (unreconstructed) fiux. The dot-dashed line 
finally represents the stabilization background 

3.4 Interactive Analysis 

After running the fitting procedure on the time history of all 
pixels making up a raster, the interactive analysis of fitting 
results is carried out, looking in detail at portions of the 
time history that were not well-fit by the automated anal- 
ysis. The details of the interactive analysis process need to 
be tuned to the quality of the specific raster under consid- 
eration. In particular, the choice of thresholds in interactive 
checks is closely related with observing parameters such as 
the exposure time but also with the varying frequency and 
severeness of cosmic ray impacts. Generally speaking, more 
careful checks could be profitably carried out on intrinsically 
higher- quality data. On the other hand, it is desirable to en- 
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Figure 3. DifTcrcnt troublesome situations in pixel time histo- 
ries: a) Recovery of stabilization background level after a fader 
b) Recovery of stabilization background level after a dipper c) 
Detection of a bright source hidden by a strong common glitch 
d) Faint source hidden by the recovery of the stabilization back- 
ground level after a dipper. Here "Background" stands for stabi- 
lization background and "Detected Flux" for reconstructed flux. 



sure the most uniform data reduction as practicable. It was 
thus decided to adopt a common trade-off between the data 
quality of all rasters making up our dataset and to apply it 
to all rasters. 

As a first step, fitting failures flagged by substantial 
data-model rms deviations (higher that 0.23 ADU/gain/s) 
or negative signals (lower than -0.6 ADU/gain/s) are 
checked. Then all sizeable signal excesses (unreconstructed 
signals higher than 0.5 ADU/gain/s) are individually in- 
spected. Whenever the need arises, a further fit extending 
to the whole pixel time history or to a smaller portion of it 
is carried out. Particularly noisy regions or very strong in- 
dividual features completely preventing data reduction are 
masked. 

The massive work of interactive analysis is carried out 
with an easy-to-use IDL widget-based Graphical User Inter- 
face, a screenshot of which is shown in Figure^] which allows 
any kind of operation that could be necessary: data visual- 
ization and browsing, glitch addition and correction, time 
history masking and re-fitting. Such software is the main 
factor allowing the increase in the v olume an d quality of the 
obtained catalogues with respect to lLari e~al.. (.200 L) . 

3.5 Mapping 

Once a satisfactory fit is obtained over the whole time his- 
tory for all pixels, one can proceed to the generation of sky 
maps and to source extraction. After masking glitches and 
other noisy parts of the time history identified during Inter- 
active Analysis, fiat-fielding is carried out and a signal esti- 
mate for each raster pointing is computed by averaging all 
readouts relative to that pointing. The result is converted 
from ADU/gain/s to mjy using ISOCAM LW3 standard 
sensitivity 1 mJy — 1.96 ADU/gain/s. The "images" thus 
computed relative to all raster pointings are then projected 
onto a 2" x 2" pixel sky map adopting nominal astrometric 
information and using a new mapping technique optimized 
to partly overcome the severe PSF undersampling due to 
the large pixel size adopted in observations. Here and in 
the following, all necessary projections are carried out using 
the projection C++ code incl uded in CIA. In so doing, 
ISOCAM severe field distortionj^kjjmyrg] 202^was taken 
into account as determined bv lAussel et al. (Il99fll|. whereas 
PSF modeling followed the prescriptions given bv lOkumura 
for stellar PSFs but adopting a spectrum of the form 
fi, = constant, i.e. a closer match to the expected galaxy 
spectrum than the Ray leigh- Jeans form used for stellar spec- 
tra. The resulting PSF is larger than the one computed for 
stars. 

Following this procedure, both unreconstructed and re- 
constructed sky maps are produced. In both cases, two an- 
cillary maps are also constructed. The sky coverage (or 
"NPIX") map contains the number of independent "images" 
added together to obtain a pixel value in the final map. On 
the other hand, the noise (or "RMS") map contains the es- 
timate of the root mean square noise of each single pixel in 
the final map, computed by scaling the overall root mean 
square noise of the map as measured at its centre according 
to the inverse square root of the noise map. 

Source detection is then performed on unreconstructed 
sky maps. First, all pixels above a conservatively low flux 
threshold (60 /iJy) are selected, then the IDL Astron- 




omy User's Library find routine (based on DAOPHOT's 
homonymous routine, particularly suited for the detection 
of point-like sources in crowded fields) is used to identify 
positive brightness perturbations around these pixels and 
return their peak flux, S/N, centroid and shape parameters 
such as roundness and sharpness. 

At this stage, further interactive checks are performed 
on all sources detected with a S/N greater than 5 to as- 
sess their reliability. More speciflcally, sources detected on 
the sky map are projected back onto the time history to 
identify all raster pointings where the source is supposed to 
contribute a significant signal. Then all these pointings are 
individually checked and, if necessary, refitted, to improve 
the fit and thus recover lost signals or remove spurious ones. 
The overall results of these further checks are an increased 
reliability and an improved astrometric and photometric ac- 
curacy. 



3.6 Mosaicing 

Up to this stage, all rasters are reduced and processed in- 
dividually. However, in order to fully exploit the limited re- 
dundancy of the observations, a technique to build a mo- 
saic out of rasters covering the same field was devised. This 
is carried out as follows. Once the reduction of all rasters 
of interest is completed according to the above procedure, 
the necessary corrections to nominal astrometry are deter- 
mined as the median offset between ELAIS source positions 
expressed with respect to nominal astrometry and the posi- 
tions of USNO A2.0 sources found in the field. This is done 
through a two-step process. First, the two catalogues are 
cross-correlated, assigning to each ELAIS sources its closest 



USNO association. The median of the positional differences 
thus determined, excluding ELAIS sources with no USNO 
association within 12 arcsec, is computed and assumed as a 
first-order correction to ELAIS source positions. The cross- 
correlation procedure is then repeated to calculate a second- 
order astrometric correction in exactly the same way, the 
only difference being that ELAIS sources with no USNO 
identification within 3 times the root mean square devia- 
tion of ELAIS-USNO association distance are also excluded 
during this second step. The sum of the two corrections is 
assumed as the raster offset with respect to nominal astrom- 
etry, and the root mean square deviation of ELAIS-USNO 
distances as the error in the offset determination. Deviations 
from nominal astrometry have long been known to be sig- 
nificant in ISO raster observations, and this is confirmed by 
our results, which are summarized in Table |3] Most total 
offsets are greater than or of the order of the pixel field of 
view size (6 arcsec), whereas associated errors have a mean 
value of 0.39 arcsec. The small errors are not only due to 
the careful data reduction, but also to the large number 
of ELAIS sources, and largely contribute to the very good 
overall astrometric accuracy quantified in Section Q 

A common mosaic grid is then built on which differ- 
ent rasters belonging to the same field are projected tak- 
ing into account astrometric offsets. Mapping, source ex- 
traction, projection of sources detected with a S/N greater 
than 5 back on time history and interactive checks are 
furtherly performed on this mosaic sky map exactly as 
on single raster sky maps. The quality of data reduction 
is thus improved through cross-checks of sources detected 
on different rasters, increasing reliability and completeness 
in repeatedly observed regions and partly overcoming the 
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otherwise severe problems at raster boundaries. The fi- 
nal S/N sky maps of the four fields are shown in Fig- 
ures 1^ El and IHl Typical noise levels are between 20 
and 30 /iJy/pixel. Such maps, together with noise and sky 
coverage maps, will soon be available in FITS format at 
http : //astro . imperial . ac .uk/ ~vaccari/elais. 



4 AUTOSIMULATION 

Even after a very careful reduction, fiux determination is a 
very delicate process owing to the interplay of different fac- 
tors, the most important being mapping effects related to 
PSF undersampling, detector's transients and other trends 
due to the limits of the adopted reduction method. Map- 
ping effects are particularly tricky, since they determine the 
way one obtain source total fluxes from corresponding peak 
fluxes measured on sky maps by source detection software. 
Clearly, the conversion factor between these two quantities 
is very much dependent on the actual source position with 
respect to the centre of the pixel, a source total flux be- 
ing nearer to its peak flux the nearer the source is to the 
center of the sky map pixel. The technique we developed, 
which we called autosimulation, was conceived to model 
and insofar as possible take into account mapping and data 
reduction effects on flux determination by first constructing 
both real and theoretical (i.e. noise-free) sky maps, and then 
use the latter ones to correct fiuxes obtained from the for- 
mer ones. Autosimulation consists in simulating sources, on 
top of noise- free maps, at each position where a source was 
actually detected (or, for simulations described in Section|Sl 
at their randomly generated positions), calculating the ratio 
between peak and total fiux for the simulated source on such 
theoretical sky maps, and converting the measured peak fiux 
of the real source to its total fiux using such ratio. Source 
simulation is carried out in a straightforward way using the 
same mathematical model the LARI method is based on, 
thus also allowing to verify its predictions. 

In order to describe how autosimulation works in some 
greater detail, a few definitions and relations which will also 
be useful to discuss the results of simulations described in 
Section 13 must be summarized: 

• fa is the measured peak flux obtained from real (i.e. 
containing glitches, noise and transients) sky maps. Its value 
therefore depends on mapping effects, transients and the 
adopted data reduction method; 

• /o is the "theoretical" peak flux obtained from simu- 
lated (i.e. containing neither glitches nor noise but taking 
into account source transients) maps. Its value depends on 
mapping effects and transients only. 

• fsr and for are analogous to fs and /o but are built from 
reconstructed maps, thus recovering the effects of transients. 
As already mentioned, this strictly holds above a certain flux 
threshold only, where the values of these quantities actually 
cease to depend on transients. 

The correction factor for mapping effects is computed as fol- 
lows. A simulated source is generated on each position where 
a real source was actually detected, assuming a total flux 
So = fs/ ifs/S)^-^ based on the measured peak flux fs and 
on the average peak flux / total flux ratio (fs/S)^-^ = 0.216 
obtained from simulations (see Section |S| and Figure [T^ . 



Then the flux estimate corrected for mapping effects is com- 
puted as S = {fs/ fo) So- Source simulation is carried out 
using the same model for charge release adopted in fltting 
the data. In so doing, however, autosimulation does not cor- 
rect for more subtle and elusive effects arising from data 
reduction, such as flux loss due to bad flts, which can only 
be assessed through simulations. This is described in Sec- 
tion |^ 

Such a procedure provides reliable flux estimates which 
are well compatible with estimates obtained through conven- 
tional aperture photometry for most sources, but is clearly 
non-optimal when dealing with extended sources. For these, 
as well as for a few very close or blended sources, aperture 
photometry usually provide a better estimate of the source 
ffux. After correcting fluxes for mapping effects, potentially 
extended sources are therefore identified through both vi- 
sual inspection of sky maps and the calculation of parame- 
ters connected to source extension. Sources fiagged by these 
criteria (totalling 67 out of 1923 making up the catalogue, 
or 3.5 %) were then individually treated, aperture photom- 
etry with a suitable aperture radius was carried out and 
the result was adopted as their most reliable fiux estimate. 
As to the overall performance of the autosimulation process 
when compared with aperture photometry, autosimulated 
fiuxes and their counterparts calculated through aperture 
photometry are compared for all sources in Section [5.31 



5 SIMULATIONS 

Due to the peculiar nature of ISO data and of the reduc- 
tion method employed, it is important to carefully test the 
performance of the latter on "ideal" data and sources. In 
particular, systematic effects on ffux estimates related to 
the data reduction method can only be probed by these 
means. Due to the strong peculiarities of our dataset, which 
is characterized by several noise features on different time 
scales, only real data can effectively be taken as represen- 
tative of instrumental behaviour. Therefore, the effects of 
artificial sources must be somehow simulated on the top of 
real pixel time histories, and data reduction must then be 
carried out exactly as done for real sources. Source confu- 
sion in the field is thus slightly increased, but this effect is 
not critical for ELAIS data. The LARI Method is based on 
a physical model of the detecto's behaviour, so can make 
straightforward predictions of the effect on the detector of 
the additional photons from simulated sources. On this ba- 
sis, an extensive set of simulations was carried out to assess 
the effects on fiux estimates and the overall performance of 
data reduction in a statistically meaningful way. The SI field 
was chosen for this purpose as the most representative, be- 
ing the largest and including both regions observed only once 
and repeatedly observed regions (hereafter, non-repeated 
and repeated regions, respectively), therefore allowing to 
assess the differences in performance warranted by higher 
redundancy. 

Simulations were carried out as follows. First a set of 
N — 600 random sky positions (excluding regions near to 
real sources, mosaic sky map boundaries and masked re- 
gions of mosaic sky map, i.e. preventing simulated sources 
to appear nearer than 15 arcsec to any of these regions) was 
generated. Then a logarithmically uniform fiux distribution 
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Table 3. Astrometric offsets of ELAIS rasters. Corrections to ISO nominal astrometry as determined through cross-correlation between 
ELAIS and USNO A2.0 source lists. 



Raster Nominal Position RA (") Dec (") Total (") 



RA (J2000) Dec (J2000) offset ± error offset ± error offset ± error 



Sl.l 


00 


30 


25, 


.40 


-42 57 00.3 


-2.06 ± 0.40 


-4.46 ± 0.38 


4.91 ± 0.55 


S1.2 


00 


31 


08, 


.20 


-43 36 14.1 


-3.24 ± 0.22 


+6.86 ± 0.29 


7.59 ± 0.36 


S1.3 


00 


31 


51, 


.90 


-44 15 27.0 


+1.57 ± 0.29 


-7.75 ± 0.33 


7.91 ± 0.44 


S1.4 


00 


33 


59, 


.40 


-42 49 03.1 


+0.23 ± 0.22 


-4.01 ± 0.27 


4.02 ± 0.35 




00 


34 


44, 


.40 


-43 28 12.0 


-3.50 ± 0.23 


+9.63 ± 0.27 


10.24 ± 0.35 


S1.5.B 


00 


34 


44, 


.40 


-43 28 12.0 


-0.52 ± 0.21 


-8.10 ± 0.26 


8.12 ± 0.33 


S1.5.C 


00 


34 


44, 


.40 


-43 28 12.0 


-3.04 ± 0.24 


+5.34 ± 0.29 


6.14 ± 0.38 


S1.6 


00 


35 


30, 


.40 


-44 07 19.8 


+0.60 ± 0.43 


-7.14 ± 0.24 


7.17 ± 0.49 


S1.7 


00 


37 


32, 


.50 


-42 40 41.2 


+1.26 ± 0.22 


-5.62 ± 0.38 


5.76 ± 0.44 


S1.8 


00 


38 


19, 


.60 


-43 19 44.5 


+0.72 ± 0.22 


-5.31 ± 0.24 


5.36 ± 0.36 


S1.9 


00 


39 


07, 


.80 


-43 58 46.6 


-2.34 ± 0.19 


+4.61 ± 0.23 


5.17 ± 0.30 


Nl 1 




1 ^ 
±0 


ni 

Ul . 


.uu 




_l_i _|_ n "^A 


Q "XA -1- n 91 


Q ^A A- C\ AC\ 




16 


13 


57, 


.10 


+54 59 35.9 


+3.27 ± 0.31 


-9.48 ± 0.22 


10.02 ± 0.38 


N1.2.B 


16 


13 


57, 


.10 


+54 59 35.9 


+4.17 ± 0.23 


-7.89 ± 0.15 


8.92 ± 0.27 


N1.3 


16 


10 


34, 


.90 


+54 11 12.7 


-1.04 ± 0.42 


-6.18 ± 0.25 


6.27 ± 0.49 


N1.4 


16 


09 


27, 


.00 


+54 49 58.7 


+4.60 ± 0.44 


-6.71 ± 0.23 


8.14 ± 0.50 


N1.5 


16 


06 


10, 


.80 


+54 01 08.0 


-8.90 ± 0.27 


+3.77 ± 0.18 


9.67 ± 0.32 


N1.6 


16 


04 


59, 


.00 


+54 39 44.3 


-8.71 ± 0.32 


+3.88 ± 0.22 


9.54 ± 0.39 


N2.1 


16 


32 


59, 


.80 


+41 13 33.2 


-4.68 ± 0.23 


+7.18 ± 0.21 


8.57 ± 0.31 


N2.2 


16 


34 


44, 


.50 


+40 38 45.0 


+2.93 ± 0.27 


-2.87 ± 0.27 


4.10 ± 0.38 


N2_3 


16 


36 


05, 


.50 


+41 33 11.8 


-5.50 ± 0.46 


+7.30 ± 0.37 


9.14 ± 0.59 


N2_4 


16 


37 


48, 


.90 


+40 58 13.1 


-5.36 ± 0.24 


+7.46 ± 0.21 


9.19 ± 0.32 


N2.5 


16 


39 


13, 


.80 


+41 52 31.6 


-5.43 ± 0.29 


+7.32 ± 0.24 


9.11 ± 0.38 


N2.6 


16 


40 


55, 


.50 


+41 17 22.7 


-2.45 ± 0.28 


+2.32 ± 0.24 


3.38 ± 0.37 


N2JI.A 


16 


35 


45, 


.00 


+41 06 00.0 


+1.29 ± 0.33 


+6.01 ± 0.21 


6.15 ± 0.39 


N2.R.B 


16 


35 


45, 


.00 


+41 06 00.0 


+5.91 ± 0.34 


+4.24 ± 0.22 


7.27 ± 0.40 


N3.3 


14 


29 


38, 


.30 


+33 24 49.6 


-1.39 ± 0.30 


-3.44 ± 0.16 


3.71 ± 0.34 


N3.5 


14 


32 


38, 


.20 


+33 11 10.3 


-0.17 ± 0.31 


-3.52 ± 0.24 


3.52 ± 0.39 



covering the range between a lower limit Si — 0.5 mjy and 
an upper limit Su = 4.0 mJy 



On = Ol \ 



n/N 



ior n = Q,N 



(5) 



was coupled with the random positions above to character- 
ize the set of simulated sources. Following the same pro- 
cedure, 100 more sources were simulated (adopting a flux 
distribution following the same analytical formula but cov- 
ering the 0.35-3.5 mJy flux range) in the repeated regions of 
the field, thus increasing the otherwise low number of sim- 
ulated sources in these regions. All 700 sources were then 
projected from their sky position back onto the pixel time 
history, simulating their effects as superposed on glitches 
and noise characterizing the real data. In so doing, all por- 
tions of time history where simulated sources were predicted 
to give a sizeable additional signal were identified. The fit- 
ting procedure was thus re-run on these regions only, and 
all steps of data reduction and flux determination were then 
carried out on simulated sources exactly as described in Sec- 
tion|21for real sources. At this stage the positions and fluxes 
of simulated sources that had been detected with a, S/N 
greater than 5 were compared with the input values to cal- 
ibrate and evaluate the performance provided by our data 
reduction method in different respects. The overall num- 
ber of simulated and detected sources is given in Table 0] 
whereeis the following Sections detail all relevant aspects of 
simulations, apart from the completeness estimates which 



Table 4. Number of Simulated and Detected Sources. Total fig- 
ures are given together with those for non-repeated and repeated 
regions. 





Non-Repeated 


Repeated 


Total 


Simulated Sources 


502 


198 


700 


Detected Sources 


230 


125 


355 



will be described by iLari et al] ^ 

extragalactic source counts. 



5.1 Flux Determination 



in order to obtain 



The autosimulation process we use for ffux determination, 
though relatively simple in principle, leads to several system- 
atic effects which need to be carefully taken into account in 
order to understand how it can provide us with the best esti- 
mate of source ffuxes. The analysis of its results on simulated 
sources will also allow us to test the goodness of our model 
and the reliability of the related signal reconstruction pro- 
cess. As described in Section |11 autosimulation involves sim- 
ulating theoretical (i.e. noise-free) maps containing all (i.e. 
both real and simulated) sources of interest, then carrying 
out signal reconstruction on both theoretical and real sky 
maps, thus determining the theoretical ratio between peak 
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Figure 9. Peak/Total Flux Ratio for theoretical sky maps. His- 
togram distribution for both Unreconstructed and Reconstructed 
maps are given. Overall distributions are shown (top panel) to- 
gether with those for non-repeated and repeated regions sepa- 
rately (bottom panel). 
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Figure 10. Real/Theoretical Peak Flux Ratio Distribution. The 
comparison between distributions for all detections and for bright 
source (> 3 mJy) detections is shown (top panel) together with 
the same comparison for detections in non-repeated and repeated 
regions (bottom panel). 



and total fluxes, and finally applying results determined on 
theoretical maps to correct real fluxes. 

Figure El shows how mapping effects affect the relation 
between peak fluxes and total fluxes for both reconstructed 
and unreconstructed theoretical sky maps. Following the 
nomenclature adopted in Section 01 the histogram distri- 
bution of fo/So and for/So ratios are plotted for all sources 
(top panel) and separately for non-repeated and repeated 
regions (bottom panel). While the two distributions for non- 
repeated and repeated regions peak at slightly different val- 
ues and thus have a (marginally) smaller width, the /01//0 
ratio is remarkably similar, 1.61 being its mean value. This 
means that, according to our model, due to the short in- 
tegration time adopted in observations with respect to the 
detector's response time constant, the measured flux, even 
in a noise-free image, will only equal 60 % of the incident 
flux. 

In order to correctly determine fluxes, however, one 
needs to make use of the employed detector's modeling and 
thus compare real fluxes with theoretical ones. By means of 
the fs/fo histogram distribution. Figure [T?71 shows how real 
peak fluxes are systematically lower than theoretical ones, 
due to limits of the data reduction method and to the pres- 
ence of noise. The width of the fs/fo distribution is quite 
large, peaking at 0.83, showing a predictable narrowing of 



the distribution at bright fluxes but on the whole a negligi- 
ble dependence on flux (Figure [Till top panel, see also Fig- 
ure ^Jl. Conversely, values for non-repeated and repeated 
regions differ in a measurable way, mean values being 0.81 
and 0.87, respectively. (Figure [Till bottom panel). The rms 
half-widths of the distributions are 0.17 and 0.12 in the two 
cases, respectively, which also provides an estimate of the 
photometric error in different regions, a quantity which will 
be assessed by different means in Section |2l 

Figure [TTI confirms the moderately small spread of the 
fs/fo distribution and also shows the good linearity dis- 
played by the fs vs. fo relation, as long as a source is de- 
tected. Provided a suitable flux correction is applied to ac- 
count for the fs/fo distribution (see Section f5.2l this allows 
to confidently use our results at all the catalogue flux levels. 

Figure fT?l shows how well the signal reconstruction pro- 
cess carried out by the fitting procedure works at different 
fluxes. At it can be clearly seen, below a certain flux, fs/fsr 
systematically falls below the fo/for, that is goodness of 
signal reconstruction breaks down at the faint end. It was 
thus decided to use only unreconstructed fluxes in flux de- 
termination, so as to provide a flux estimate that would be 
reliable at all flux levels. 

Finally, and most importantly. Figure ITsl shows how to- 
tal injected fluxes and peak real fluxes are related to each 




Figure 11. Real vs. Theoretical Peak Flux. The fs = 0.83 /o 
line, corresponding to the 0.83 mean value of the fs/ fo ratio for 
detected sources, is also plotted. 



Figure 14. Flux Correction applied to Autosimulated Fluxes. 
The inverse of the correction factor, or Qmed^ is plotted against 
S/N for non-repeated and repeated regions. 
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Figure 12. Signal reconstruction as function of real peak flux. 
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Figure 13. Total vs. Peak Flux. The S = 4.64 /s line, corre- 
sponding to the 4.64 mean value of the S/fs ratio for detected 
sources, is also plotted. 



other. The relation mimics the one between fs and fo al- 
ready shown in Figure 1111 a linear relation with a reason- 
ably small spread (the mean value of the S/fa ratio being 
4.64), warranting a reliable determination of total fluxes on 
the basis of measured peak fluxes. 



5.2 Flux Correction 

Systematic effects on flux estimates described in previous 
Section were modeled and taken into account using results 
from simulations and following Gruppio ni ct al . ( 2002). Ac- 
cording to this approach, the correction to be applied to 
flux estimates obtained through autosimulation is derived 
from the the so called g function, which describes the S/N- 
dependent distribution of the fs/ fo ratio. The g function is 
obtained as the combination of the intrinsic (i.e. high-5'/A'^) 
g function, or go function, and a term due to noise. First, 
the go function is derived from the distribution of the fs/ fo 
ratio obtained for bright (S > 3 mjy) simulated sources (see 
Figure by modeling and correcting for the small degree 
of incompleteness to be expected at such fluxes. Then con- 
volution with a variable noise term is carried out to obtain 
the overall S/A'^-dependent g function. Finally, for each S/N 
value, the median fs/ fo ratio for detectable (i.e. fs > 5a) 
sources, or Qmed, is computed. Autosimulated fluxes are then 
corrected by a factor 1/qmed- This process was carried out 
separately for non-repeated and repeated regions, and the 
results are plotted in Figure IT^ The two Qmed curves asymp- 
totically tend to the average fs/ fo values 0.81 and 0.87 de- 
termined for simulated sources in non-repeated and repeated 
regions, then virtually coincide and soar below S/N ~ 7, 
where qmed ~ 1. 

It must be emphasized how all of the above applies to 
"relative" fluxes determined with respect to ISO LW3 stan- 
dard sensitivity. In Section|5|the absolute flux calibration of 
our catalogue will be discussed on the basis of the compari- 
son between measured and predicted MIR stellar fluxes and 
with respect to IRAS standard calibration. 
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Figure 15. Autosimulated Fluxes vs. Aperture Fluxes 



5.3 Autosimulation vs. Aperture Photometry 

The performance of the autosimulation technique for flux 
determination described in Section |1] can also be assessed 
comparing the total fluxes obtained with this procedure with 
those obtained by means of ordinary aperture photometry. 

The relation between these two quantities shows a good 
linearity for most sources, as it is shown for simulated 
sources in Figure [T^ Plotted aperture fluxes were computed 
adopting an aperture of 6 arcsec radius and correcting both 
for the transients and for the 40 % of the instrumental PSF 
falling beyond this aperture. Such an aperture was chosen as 
the most reliable trade-off allowing to reliably include most 
of the source flux and lest of the background. 

Only a few sources depart substantially from the 1:1 re- 
lation, with 94 % of sources with fluxes in accordance within 
20% and an overall rms deviation of 12%. On the whole, 
our flux determination procedure is therefore consistent with 
conventional aperture photometry, provided a proper correc- 
tion for PSF effects is applied. 

5.4 Catalogue Flux Distribution 

Flux distribution of catalogue sources is illustrated in Fig- 
ure 1161 with 50 % of the sources having fluxes greater than 
1.65 mJy, 76% greater than 1.2 mJy and 89% greater than 
1 mJy. As it can be clearly seen from the two separately 
plotted histograms, sources in non-repeated regions largely 
prevail in number down to about 1.5 mJy, at which flux their 
number per flux bin drops sharply. Conversely, the number 
of sources per flux bin in repeated regions continue to in- 
crease, if slowly, down to 1 mJy, where they already amount 
to about 50% of the total number of sources. While de- 
tailed completeness estimates will be described bv lLari et alJ 
i2004fl . this plot seems to indicate that the catalogue is es- 
sentially complete up to ~ 1.5 mJy, and slightly fainter than 
that in repeated regions. 



6 OPTICAL IDENTIFICATIONS 

Identification of 15 fim sources was carried out on both 
archival and deep imaging obtained for this purpose at op- 
tical and near-infrared wavelengths. This provided optimal 
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Figure 16. Histogram Flux Distribution of ELAIS Sources. All 
regions, non-repeated and repeated regions are plotted. 



stellar discrimination as well as optical and near-infrared 
flux measurements over the whole large flux range probed 
by our catalogue, thus ultimately leading to a very reliable 
15 pm photometric calibration (see Section |UJ. 

The SI field was surveyed in R down to 7? ~ 23.0 by 
iLa Franca et al.l (|20o3) using the 1.5 Danish/ESO telescope. 
The Nl and N2 fields were surveyed in U, g' , r', i' and Z 
down to 23.4, 24.9, 24.0, 2 3.2 and 21.9 respectiv ely, as part 
of the Wide Field Survey jMcMahon et al.ll200ll) . using the 
Wide Field Camera at the Isaac Newton Telescope. Identi- 
fications of Nl and N2 15 ^, m sources based on maximum 
likelihood were obtained by llConzalez-Solares et al.ll2004f) . 
Deep near-infrared imagin g of areas around sele cted 15 pm 
sources was carried out by llVaisanen et al .1200211 . Furtherly, 
USNO A2.0, CSC 2.2, Tycho-2, 2MASS AU Sky Data Re- 
lease and APM catalogues were cross-correlated with 15 pm 
sources so as to provide a list of tentative identifications. 15 
^m contours were superimposed on DSS and 2MASS images 
as well as on new observational material, yielding diagnostic 
finding charts for extended, blended and disturbed sources. 

Generally speaking, the identification beised on deep 
optical imaging and automated source extraction and clas- 
sification carried out with SExtractor was usually chosen, 
but archival material proved essential in dealing with bright 
sources, which appeared as saturated in deep observations, 
and in providing measurements for sources in the N3 or in 
small portions of the other three fields, where deep optical 
imaging had not been obtained. In particular, a positional 
difference of 6 arcsec, corresponding to 3 times the maximum 
estimated astrometric error of the 15 pim data reduction pro- 
cess (see SectionlSJ, was adopted as a cut-off ensuring a safe 
identification. 

The magnitude distribution of the sources identified as 
optical counterparts of 15 /xm sources through this process 
is plotted in histogram form in Figure lTTI showing a bimodal 
distribution whose two peaks are to be associated with stars 
and galaxies and sharply cutting-off at ii ~ 23. The number 
of stars and total sources identified according to this crite- 
rion, is separately reported for the four fields in Table |S| 
About 90 % of 15 ^m sources are assigned a robust optical 
association in Nl and N2 fields, while due to the shallower 
depth this fraction drops to about 80% in the SI and N3 
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Figure 17. Magnitude distribution of optical IDs of 15 
sources. The two peaks at r ~ 12 and r ~ 18 are associated 
with Stellar and Extragalactic IDs, respectively. 



Figure 18. 15 fira vs. R Flux-Flux Diagram of Identified Sources. 
Stellar and Extragalactic IDs are indicated. 



Table 5. Some basic properties of the catalogue divided into 
different fields. 



Field 


Area 


Sources 


Stars 


IDs 




[deg2] 


# 


# 


% 


# 


% 


SI 


4.17 


736 


145 


19.7 


584 


79.3 


Nl 


2.84 


490 


121 


24.7 


441 


90.0 


N2 


2.84 


566 


126 


22.3 


493 


87.1 


N3 


1.00 


131 


29 


22.1 


103 


78.6 


Total 


10.85 


1923 


421 


21.9 


1621 


84.3 



fields The average fraction of identified sources amount to 
84.3 %. Conversely, the statistics of stellar identifications is 
remarkably uniform between the 4 fields, with a 22 % of stars 
and small field-to-field differences. 

Figures [TSl and lT^ show how, at least down to the 15 /im 
flux limit probed by ELAIS, it is furtherly possible to effec- 
tively discriminate between stars and galaxies solely on the 
basis of the comparison between optical and near-infrared 
magnitudes with 15 jim. fluxes, thus demonstrating the ro- 
bustness of the identification process. 

The optical identification process also allowed us to as- 
sess the astrometric accuracy achieved in 15 /im data re- 
duction independently of simulations, as it is described in 
Section □ 



7 ASTROMETRIC ACCURACY 

Long-wavelength observations are frequently hampered by 
large astrometric errors, leading to severe difficulties in 
multi- wavelength identifications and thus physical studies of 
detected sources. In our case, astrometric errors can be con- 
sidered as the combination of three error sources associated 
with the detector spatial sampling Os , the reduction method 
ffr and the instrumental pointing accuracy ap, respectively. 

The combination of the CTs and ar terms, which we will 
hereafter indicate as a^+r, can be evaluated from simula- 
tions, comparing the injected positions of simulated sources 
with the corresponding detected positions, whereas the ap 




10 12 



R - [15] 



Figure 19. R — K vs. R — [15] Color-Color Diagram of Identified 
Sources. Stellar and Extragalactic IDs are indicated. The [15] 
magnitude scale is defined in Equation 151 



term is given by the error on the rasters' astrometric offset, 
as derived in Section 13.61 The total astrometric error will 
then be given by 



a tot 



(6) 



The overall distribution of the differences between injected 
and detected positions of simulated sources is shown in his- 
togram form in Figure 1201 for both non-repeated and re- 
peated regions. Differences in RA and Dec are distributed 
in a remarkably similar way, which will allow us to use them 
together in the following statistical analysis, so as to in- 
crease our sample. About 75 % simulated sources are de- 
tected within 1 arcsec of their injected positions and about 
96 % within 2 arcsec. Dependence of the Oa+r on SjN can be 
evaluated binning simulated sources according to their SjN 
and computing the median positional difference of each bin. 
Thus one obtains the Os-vr term as function of SjN . The ap 
term is obtained as the median value of astrometric offset 
errors per component reported in Tabled that is 0.39 arcsec 
per component irrespectively of 5'/A^. The total astrometric 
error per component is then computed through Equation |S| 
The results are illustrated in Figure l^ showing total astro- 
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Figure 20. Sampling and Reduction Astrometric Error. Posi- 
tional differences between injected and detected simulated sources 
in RA and Dec. 
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Figure 21. Astrometric Accuracy. Total astrometric errors per 
component as function of S/N in non-repeated and repeated re- 
gions. Polynomial fitting curves used for the evaluation of cata- 
logue astrometric error are also indicated. 



metric error begin dominating the astrometric error budget, 
and the astrometric error curve starts showing irregularities, 
so that linear fitting cannot be reliably assumed to describe 
the actual astrometric error. Accordingly, values obtained 
from simple extrapolation of the straight line obtained at 
S/N < 20 to higher 5'/A'^ were truncated when they were 
deemed too optimistic, i.e. at the 0.8 and 0.7 arcsec values 
for non-repeated and repeated regions, corresponding to the 
extrapolated values for S/N ~ 25. 



metric errors for sources detected both in non-repeated and 
repeated regions as function of S/N for S/N < 30. The as- 
trometric accuracy turns out to be very good, errors being as 
low as 0.7 and 0.8 arcsec per component at the bright end in 
non-repeated and repeated regions, respectively. Catalogue 
astrometric errors can be estimated by fitting polynomials 
to the log atot vs. log S/N curves for non-repeated and re- 
peated regions derived from simulations. At the faint end, 
however, the errors computed in repeated regions unrealis- 
tically flatten below S/N ~ 10, due to the poor statistics 
following from the low number of simulated sources that are 
detected at low S/N levels. For this reason, we cross-checked 
the astrometric accuracy evaluated from simulations against 
an independent estimate based on optical identifications of 
catalogue sources (see Section [SJ. The same binning process 
described above was carried out on positional differences 
between catalogue sources and their optical counterparts, 
yielding the results plotted in Figure 

HI] The thus-derived 
(Jtot vs. S/N relation closely follows a straight line in the 
log-log plane at 5 < S/N < 20, and closely resembles the 
one derived from simulations at 10 < S/N < 20, confirming 
the overall robustness of the simulation process. It was thus 
decided to estimate catalogue astrometric errors by fitting 
a straight line to the log atot vs. log S/N curves for non- 
repeated and repeated regions derived from optical identifi- 
cations. At S/N > 20, other effects such as the widening of 
the instrumental PSF and the increase of the optical astro- 



8 PHOTOMETRIC ACCURACY 

Errors in flux determination using our method can mainly be 
attributed to two effects, namely the autosimulation process 
and the noise present on the sky maps. The first contribu- 
tion can be estimated by computing the width of the /s//o 
distribution shown in Figure [THl for high S/N sources only, 
so as to evaluate the effects of the autosimulation process on 
relatively noise-free maps. At lower S/N, photometric errors 
arising from noise on the sky maps will combine with those 
arising from autosimulation. The overall error is thus given 
by 




where A{fs/fo) is the width of the fs/fo distribution as 
measured for high S/N sources for non-repeated regions and 
repeated regions separately, and amap is the noise as mea- 
sured at each source position on the sky map. The first term 
is a constant, about 0.15 and 0.11 for non-repeated and re- 
peated regions, respectively. This term dominates the photo- 
metric error budget at high S/N, whereas the importance of 
the second term increases when S/N decreases. Photometric 
errors estimated using Equation |7| are shown in Figure 1^ 
The values for S/N ~ 10 are 0.18 / 0.15 for non-repeated 
/ repeated regions, while they increase to 0.25 / 0.23 for 
S/N ~ 5. 
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Photometric Calibration 




Figure 22. Relative photometric error as function of S/N as 
given by Equation |7| for non-repeated and repeated regions. 



9 PHOTOMETRIC CALIBRATION 

The standard sensitivity of ISO LW3 channel is of 1.96 
ADU/Gain/s/mJy, which can be used to cahbrate our cata- 
logue fluxes. However, given the large sample at our disposal, 
the accuracy of such photometric calibration can be tested 
against the independently determined IRAS standard pho- 
tometric calibra tion by studyi ng detected sources with stel- 
lar counterparts . I Aus'sell i2004) performed a detailed study of 
mid-infrared emission from stars, using large samples drawn 
from IRAS Faint Sour ce Catalog with counterparts in the 
2MASS and Tycho-2 llHog et al-lboOOt) catalogues. In the 
two cases, it is found that the J — K and B — V colors of 
stars are extremely well correlated with the K — [12] and 
B — [12] colors, respectively, where [12] is a magnitude scale 
constru cted from IRAS 12 u m flux, following the prescrip- 
tions of lOmont et alJ 1^93)- This relation allows to accu- 
rately predict the IRAS 12 /xm flux of a star, provided that 
its J and K (or B and V) are known and that J — K (or 
B — V) a re within certain limits. Stellar atmosphere mod- 
els IILeieune et al. 1998) show that such color criteria select 
stellar spectral types for which the ratio between the 15 fim 
flux and the 12 fim flux is essentially constant. Thus, if one 
has access to a substantial number of optical and/or near- 
infrared counterparts to 15 /im sources, it is possible to use 
the correlations bv .Aussel ( 2004) to predict the 15 /xm fluxes 
of ISO-detected stars and compare them to the measured 
values. Such a comparison then allows to investigate possible 
systematic differences between the two sets of fluxes, which 
should be attributable to discrepancies between the inde- 
pendently established IRAS and ISO calibrations. Given the 
tighter nature of the near-infrared/mid- infrared correlation 
over a wide range of fluxes with respect to the optical/mid- 
infrared one, it was decided to use the former, that is ex- 
pressed by 



K ~ [15] = 0.044 + 0.098 {J - K) , 
where [15] is a magnitude scale deflned as 
[15] = 3.202 - 2.5 log (S15 [mjy] ) , 



(8) 



(9) 



to compute predicted fluxes. 2MASS All Sky Release fluxes 
were used for the majority of the sources together with fluxes 
determined for fainter sources as part of the near-infrared 



Figure 23. Photometric calibration. Measured vs. Predicted 15 
/im Stellar Fluxes. The 1.0974 correction factor to measured 
fluxes is here already taken into account, and the one-to-one re- 
lation is also plotted. 




Meosured Flux / Predicted Flux 

Figure 24. Photometric calibration. Histogram Distribution of 
Measured / Predicted 15 ^tm Stellar Flux Ratio. The 1.0974 cor- 
rection factor to measured fluxes is here already taken into ac- 
count. 



follow -up program to ELAIS observations aisan en et alJ 
l2002h . In Figure 1221 measured fluxes are plotted against pre- 
dicted ones for all 408 sources with reliable stellar identifica- 
tion and near-infrared magnitudes, showing the remarkable 
linearity and the small spread of the relation over a wide 
range of fluxes. Figure 1^ shows the histogram distribution 
of the measured/predicted flux ratio. Reducing the sample 
to the 300 sources with a S/N > 10 yields an average pre- 
dicted/measured ratio of 1.0974 with a standard deviation 
of 0.0121. The difference between measured and predicted 
fluxes is significant at ~ 8 sigma, suggesting the presence of 
systematic effects in either the IRAS or ISO calibration pro- 
cess. Disentangling the effects leading to such a discrepancy 
would involve the reduction of a wide set of IRAS-detected 
sources observed by ISO using our method. Given the sub- 
stantial amount of work that would be needed for this pur- 
pose, it was decided to simply correct our catalogue fluxes 
by a constant 1.0974 factor to put them on the IRAS scale. 
However, while this choice was taken for the sake of compati- 
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bility of our results with studies adopting the IRAS flux scale 
(and particularly IRAS-based source counts and luminosity 
functions), this is not to indicate that IRAS calibration is 
more secure that ISO's. In fact, results from data reduction 
of deep ISOCA M surveys in the Lockman Hole using the 
LARI Method (|Fadda et al.ll20n4) favour the ISO calibra- 
tion. On the basis of stellar atmosphere models applied to a 
smaller number of stars and deep multi-band imaging, the 
authors find that model predictions agree with 15 ^m fluxes 
determined following ISO calibration. 



10 THE CATALOGUE 

The ELAIS 15 /im Final Analysis Catalogue (Version 1.0) 
contains 1923 sources detected with a S/N greater than 5 
in the (RA-ordered) four fields SI, N3, Nl and N2, totalling 
an area of 10.18 deg'^. For each entry, the catalogue reports 
astrometric and photometric information, optical identifi- 
cation and a number of ancillary flags. The entries are as 
detailed in the following; 

• Name : lAU source name constructed from Right As- 
cension and Declination; 

• RA (J2000) : Right Ascension at epoch J2000 in both 
decimal and sexagesimal angular units; 

• Dec (J2000) : Declination at epoch J2000 in both dec- 
imal and sexagesimal angular units; 

• Total Flux : source total flux obtained from autosim- 
ulation or aperture photometry (see "Aperture" below), ex- 
pressed in mjy; 

• Peak Flux : source peak flux measured on unrecon- 
structed maps, expressed in ^Jy/pixel; 

• S/N : signal-to-noise ratio measured on unrecon- 
structed maps 

• Astrometric Error : Astrometric error as determined 
from Equation |H] 

• Photometric Error : Photometric error as deter- 
mined from Equation [7| 

• Aperture : aperture photometry flag, for sources 
whose flux was estimated through autosimulation and 1 for 
sources on which aperture photometry was preferred (see 
Section ; 

• Optical ID : optical identification as described in Sec- 
tion 1^ S stands for star, G for galaxy and * for unidenti- 
fied source. A further D flags a particularly dubious (non- 
)identification. 

A sample of the catalogue, which will shortly become avail- 
able at http://astro.imperial.ac.uk/~vaccari/elais, 
is shown in Table 15] Some of the data on the sources are 
also contained within the ELA IS band-merged catalogue by 
iRowan-Robinson et alJ i2003l) . 



more robust and, most importantly, quicker and easier to ap- 
ply to large datasets. Its application, in the new form, to the 
four fields composing the dataset (inclu ding a r e-rcduction 
of SI observations already presented in iLari et a l. (2001)) 
has produced a catalogue of 1923 sources spanning the 0.5 
- 100 mJy range, detected with a S/N greater than 5 over 
a total area of 10.85 deg^. Optical identification of 15 fim 
sources has been carried out on heterogeneous optical and 
near-infrared imaging material, allowing to determine a ro- 
bust association for about 85 % of the sources and identify 
22 % of them as bona fide stars, furtherly demonstrating the 
reliability of our data reduction process. 

The evaluation of the catalogue's quality has been 
carried out through both accurate simulations and multi- 
wavelength identification. The astrometric accuracy is of or- 
der 1 arcsec in both RA and Dec for S/N > 10, while it 
increases up to about 2 arcsec in both RA and Dec for 
S/N ~ 5, and somewhat better for sources detected in 
higher-redundancy sky regions. The photometric accuracy 
is estimated to be below 25 % over the whole range of fluxes 
and redundancy levels probed by our catalogue, and better 
than 15 % for S/N > 10 sources. 

The comparison of measured stellar fluxes with fluxes 
estimated on the basis of stellar atmospehere models cali- 
brated on IRAS data and on near-infrared photometry al- 
lowed to achieve an IRAS /ISO relative photometric calibra- 
tion. An IRAS/ISO relative calibration factor of 1.0974 ± 
0.0121 was determined, shedding doubts on the goodness of 
the two independently determined calibrations at the 10 % 
level. For lack of a simple way to identify error sources in 
IRAS and/or ISO calibratioon process, it was decided to put 
our catalogue on the more c ommonly used IR AS flux scale. 

In a forthcoming paper jLari et alJl2004) completeness 
estimates and extragalactic source counts from this cata- 
logue will be presented, covering the crucial flux range 0.5 - 
100 mJy between ISOCAM 15 ^m Deep Surveys and IRAS 
All Sky Survey. 
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11 CONCLUSIONS 

A technique for ISO-CAM/PHOT data reduction, the LARI 
method, was variously refined and applied to ELAIS 15 //m 
observations. The mathematical model for the de tector's be- 
haviou r is the same as originally presented in iLari et all 
i200ll) . but thanks to various improvements, and particu- 
larly to a new Graphical User Interface, the method is now 
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Table 6. A sample portion of ELAIS 15 /^m Final Analysis Catalogue. Source lAU Name, Right Ascension and Declination at J2000 
epoch in sexagesimal formats, Total and Peak Flux, Signal- To-Noise Ratio, Astrometric and Photometric Errors, Aperture Flag. Note 
that a few columns described in text were dropped for the sake of brevity. 
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